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1. Introduction 

Propagators — that is Feynman Integrals depending on only one external momenta — ap- 
peared in QFT from its very beginning, for instance, in the description of the vacuum polarization 
in QED and since then constitute an important set of Feynman integrals to deal with. 

Massless propagators are an indispensable tool in Renormalization Group (RG) calculations 
within the framework of Dimensional Regularization [jl|] and the so-called Minimal Subtractions 
schemes 

Massless propagators appear in many important physical applications. The total cross-section 
of e^e^ annihilation into hadrons, the Higgs decay rate into hadrons, the semihadronic decay rate 
of the T lepton, the running of the fine structure coupling are all computable in the high energy 
Unlit in terms of massless propagators. 

The strong coupling constant is one of the three fundamental gauge couplings constants of 
the Standard Model of particle physics. Its precise determination is one of the most important aims 
of particle physics. One of the most precise and theoretically safe determination of is based on 
measurements of the cross section for electron-positron annihilation into hadrons. These have been 
performed in the low-energy region between 2 GeV and 10 GeV and, in particular, at and around 
the Z resonance at 91.2 GeV. Conceptually closely related is the measurement of the semileptonic 
decay rate of the T-lepton, leading to a determination of at a scale below 2 GeV. 

From the theoretical side, in the framework of perturbative QCD, these rates and cross sec- 
tions are evaluated as inclusive rates into massless quarks and gluons [^, (Power suppressed 
mass effects are well under control for g+e^ -annihilation, both at low energies and around the Z 
resonance, and for T decays [||, |6|, ^, ^, ^ |T^, and the same appUes to mixed QCD and electroweak 
corrections [jlT], p^]). 

The ratio R{s) = o{e^e^ hadrons) /a (e+e^ l^^l^- ) is expressed through the absorptive 
part of the correlator 

n^v(<?) = // dx^''^"(o|r[ 77^77(0) ]|o) = (-g^v<?'+^M^v)n(-<?2), d-D 

with the hadronic EM current j^™ = Y^fQfWfYnYf^ ^^'^ Qf being the EM charge of the quark /. 
The optical theorem relates the inclusive cross-section and thus the function R{s) to the disconti- 
nuity of n in the complex plane 

R{s) = 12 71 imn{-s-i5). (1.2) 

The renormalization mode of the polarization operator Tl{Q^) reads (see, e.g. [^) 

l+e^n{Q^/n\a,)=zf + e^no{Q\a^),, (1.3) 

where e = \faATl and a is the (renormalazied) fine structure constant. Eq. ( |1.3[ ) can be naturally 
deduced from the connection between 11(2^) and the photon propagator D^y{q) 
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It is also convenient to introduce the Adier function as 



D{Q^) = -nn'Q^^Yl^Q^) = / (1.4) 



We define the perturbative expansions 

oo oo 

D{Q^) = £ diaiiO"), R{s) = £ m[{s), (1.5) 

!=0 (=0 

where Us = as/ n and we have set the normalization scale = or o jj.^ = s for the Euclidian and 
Minkowskian functions respectively. The results for generic values of n can be easily recovered 
with standard RG techniques. 

For the vector correlator the terms of order and have been evaluated nearly thirty and 
about fifteen years ago |13, 14], respectively. The corrections are conveniently classified ac- 



cording to their power of rif, with Uf denoting the number of light quarks. The a^«y- term is part 
of the "renormalon chain", the evaluation of the next term, of order a'^Uj, was a test case for the 
techniques used extensively in this paper and, furthermore, led to useful insights into the structure 
of the perturbative series already [|^. The rest: two remaining most difficult terms of orders a^n^- 
and a^rpj- of equivalent complexity' are not yet known. The results of their (partial) evaluation will 
be presented in the talk. 

2. Calculation of <i4 at = 3 

The complete five-loop calculation requires the evaluation of 2671 Abelian quenched plus 
about seventeen thousand non-abelian and non-quenched diagrams (see Fig 1). 




Figure 1: (a) (a),(b) Examples of five -loop quenched Abelian diagrams of order a'^n'j- . (c) A five -loop 
non-abelian diagram, (d) A five -loop non-quenched abelian diagram of order a'^n f 

Using "infrared rearrangement" [|l^], the R* operation [17] and the prescriptions formulated in 



[ |18[ ] to algorithmically resolve the necessary combinatorics, it is possible to express the absorptive 



part of the five-loop diagrams in terms of four-loop massless propagator integrals ("p-integrals"). 



Using then a representation for Feynman integrals proposed in [ ]19| , [20| ], these p-integrals can be 
reduced to a sum of 28 master integrals, with coefficients which are rational functions of the space- 
time dimension D. These coefficients were, in a first step, evaluated in the large-D limit, and. 



In the sense that they both contain all possible topologies. 
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after evaluating sufficiently many terms of the 1/D expansion, fully reconstructed |21]. This direct 
and largely automatic procedure required enormous computing resources and was performed on a 
32+8-node SGI parallel computer with SMP architecture and on the HP XC4000 supercomputer of 



the federal state Baden- Wiirttemberg using a parallel version [ |22[ ] of FORM [|23|]. 

The construction of the large-D limit requires in general huge storage resources, which nat- 
urally constrains the structure of the input p-integrals: they should better not contain any extra 
parameters Uke color coefficients, ?i/-the number of Ught quark flavours contributing to internal 
fermion loops-, and so on. As a result we are forced to use a "slice" approach: that is to set all 
color coefficients to their numerical values and to fix to some integer. Thus, in order to compute 
the remaining two not-yet-known qf contributions to R{s) one should compute two slices: rif = 
and rif = hq. Here no stands for any non-zero integer which could be chosen at will. Since the 
evaluation of every single slice is a problem by itself, we decided to start from hq = 3 as the result 
has important physical applications for the analysis of the QCD corrections to hadronic T-decays 
(see, e.g. 

Our result reads (we suppress the trivial factor 3 L/2/ throughout^) 

78631453 1704247 ^ 4185 ^2 34165 ^ 1995 ^ 

^^^"^=^^ = ^0736 43^^^ + ^^^+^^^-^^' ^^-'^ 

^ 49.0757. (2.2) 

The corresponding expression for R{s) is: 

l+a,+ 1.6398a,2 + 6.3710a? -106.8798a^ (2.3) 

Since it will presumably take a long time until the next term of the perturbative series will 
be evaluated, it is of interest to investigate the predictive power of various optimization schemes 



empirically. Using the principles of "Fastest Apparent Convergence" (FAC) [24] or of "Minimal 



Sensitivity" (PMS) [pq], which happen to coincide in this order, one gets [26, 27] 

df^'^{nf = 3) =27 ±16, 
with the central value of differing significantly from the exact result 

^/f^='(„^. = 3)=49.08. (2.4) 



However, within the error estimates [|270, predicted and exact values are in agreement. The picture 
changes, once these estimates are used to predict the coefficient r4. The prediction for the final 
result is significantly closer (in a relative sense) to the result of the exact calculation: 

rf^'^inf = 3) = -129 ± 16, rl^^'\nf = 3) = -106.88. (2.5) 

This is in striking contrast to the case of the scalar correlator, where the predictions for the dy- 
namical terms work well, but, as a consequence of the strong cancellations between dynamical and 



kinematical terms fail completely in the Minkowskian region [|28[]. 



In this paper we present the results for the so-called "non-singlet" diagrams, where one and the same closed 
quark line is connected to the external currents. These are sufficient for a complete description of T-decays. For e+e^ 
annihilation they correspond to the dominant terms proportional f Of - The singlet contributions proportional (L/ 2/)^ 
arise for the first time in ffi^a^). They are known to be small, and will be evaluated at a later point. 



4 



Massless propagators: 



K. G. Chetyrkin 



3. Colour structure of 

In general, it would be of great interest to know the contribution for a generic color group. 
The structure of the result for, say, Adler function, may be predicted from [^] : 

t/4 = ^ d4,iC4j, (3.1) 

!=1,12 

2 

with the coefficients ^4 , being polynomials in In(^) and the following colour factors read 



C4 — {duCp , dRCpCA , dRCpC^ , dRCfC^ , dRCpTpnf , dRCpCATprif , 
dRCpClTpnf, dRCjT^nj , dRCpCAT^n} , dRCpT^n) , df^^df"^ , nfdf'^df'^}. 



(3.2) 



Here [T"T'']ij = Cfdij and f"<^'^f'"^d _ (j^§ab quadratic Casimir operators of the fundamental 

and the adjoint representation of the Lie algebra and tr{T''T^) = Tpd"^ is the trace normalization of 
the fundamental representation. dR is the dimension of the fermion representation (i.e. the number 
of quark colours) and rif is the number of quark flavors. The exact definitions for the higher order 
group invariants df"^df"^ and df'^'df''' are given in [^]. For QCD with SU(3): 

Cf = 4/3 , Ca = 3 , = 1/2 ,dR = 3, df"^"^ d'}^'^'^ = — , df"^'^df"^'^ = — . 
We observe that to get the full colour structure of J4 we should compute as many as eight extra 



slices (as four slices are already available: ^410 pOO, J4.9 and J4 8 []15[], as well as SU{?>) with 
rif = 3). Thus, complete evaluation of J4 for generic colour group is a matter of future. 

However, not all slices are equally difficult. For the case of the abelian U(l) group (QED!) the 
colour coefficients have especially simple form: 

Cp = l,CA=0,Tp = l,dR = l, df'^'df"^ = 0,df'''df"' = 1. 

In addition the corresponding diagrams are somewhat simpler due to the absence of the three- and 
four-gluon couplings as well as any vertexes with ghost lines. An even simpler slice corresponds 
to case of quenched QED (qQED): 

Cp = l,CA = 0,TF = l,dR = l,nf = 0,df'''df''^ = 0,df"^df"^ = 0. 

Thus, in general the quenched QED contribution to dn is equal to the coefficient in front of the 
dR {CfY structure in the colour decomposition of <i„. 



Up to order a the Adler function for qQED reads [14| 



3 fiQ 

DiQ™(a,e) = l + 3A-^A2-^A^ (3.3) 

where A = A{}i) = Note that Di^™(a,/i/2) does not depend on its second argument at 

all, which means that it is an essentially finite, scheme-independent quantity. It is instructive to 



compare eq. (3.3 ) to the j8-function of qQED, also known to the same order from [31, 32, 33|: 
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One observes the validity of the relation 



which can be easily proven in all orders of perturbation theory starting from ( |1.3[ ). Thus, the 
evaluation of the five-loop contribution to ^^'^^^ provides with us another slice, namely with d^^\. 

On the other hand, such a calculation is of importance by itself. Indeed, the j8- function of 
qQED possesses a number of remarkable and intriguing features. 

• It is gauge and scheme independent in all orders. 

• Its coefficients are simple rational numbers at three and four loops^. 

• If ^°'^^^{(Xq) = then a = Oq leads to a self-consistent finite solution of (massless) QED 
|33 



A detailed knot-theoretic explanation of the rationality property at three loop level was given 



in [36, 37]. At four loop the problem was thoroughly investigated with the help of a dimensionally 



background-field method in [|^, |3^. Some hope that the rationality property is not accidental 



but also holds in higher orders has been expressed in []38|]. Unfortunately, no clear unambiguous 
prediction as for the structure of higher (five loops and beyond) orders has ever been made. 
We have computed di\^\. Our five loop result for the j8 function reads: 

^qQED^4^^4^2_2^3_46^4^|^l^ ^ 128C3)a5. (3.6) 

Thus our calculation of dn^\ has solved the problem: the rationality ceases to exist starting from 
fifth loop. 

Acknowledgments. We thank David Broadhurst and Dirk Kreimer for useful discussions. 
This work was supported by the Deutsche Forschungsgemeinschaft in the Sonderforschungsbere- 
ich/Transregio SFB/TR-9 "Computational Particle Physics", by INTAS (grant 03-51-4007) and 
by RFBR (grant 05-02-17645). The computer calculations were partially performed on the HP 
XC4000 super computer of the federal state Baden- Wiirttemberg at the High Performance Com- 
puting Center Stuttgart (HLRS) under the grant "ParFORM". 

Note added. Recently we have finished the calculation of d^ for a generic value of nf. The 



reader is referred to [39] for details and for the discussion of phenomenological applications of the 
result. 
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